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Abstract 



We derive an expression for the effective second-harmonic coefficient of a 
dilute suspension of coated spherical particles. It is assumed that the coat- 
ing material, but not the core or the host, has a nonlinear susceptibility for 
second-harmonic generation (SHG). The resulting compact expression shows 
the various factors affecting the effective SHG coefficient. The effective SHG 
per unit volume of nonlinear coating material is found to be greatly enhanced 
at certain frequencies, corresponding to the surface plasmon resonance of the 
coated particles. Similar expression is also derived for a dilute suspension of 
coated discs. For coating materials with third-harmonic (THG) coefficient, 
results for the effective THG coefficients are given for the cases of coated 
particles and coated discs. 



I. INTRODUCTION 



With the advancements in nanotechnology, it has become possible to fabricate nanopar- 
ticles of various kinds with specific geometries. For years, random or ordered composite 
materials consisting of two or more materials with different physical properties have been 
studied intensively, with the aim of tuning the effective physical properties of the compos- 
ite by properly choosing the constituent materials, the structure of the composite, and/or 
the volume fraction of each of the constituents. Ordered composites, for example, have led 
to the development of semiconductor heterostructures and devices, and photonic band-gap 
materials. Composites consisting of small particles of one material randomly embedded in 
a host medium also show interesting behavior [1] . The percolation effect that occurs as the 
volume fraction of one component increases leads to a qualitative change in the physical re- 
sponse of a random composite. In the dilute concentration limit, various physical responses 
are affected by local field effects due to the inhomogeneous nature of the system. 

Of particular interest is the optical response of nanoparticles and nanoparticles with 
specially designed geometry [2-5]. For example, using a nanoshell geometry on particles, 
it has been found that the effective Raman scattering can be enhanced by as much as 10 6 
[6]. The surface plasmon resonance is found to be shifted in nanoparticles coated with gold 
[7]. The optical properties of nanoparticles may also be affected by quantum size effects [8]. 
In general, it is expected that the local field effects play a significant role in the nonlinear 
response in composites consisting of small particles as the response depends on high powers 
of the field [9]. Greatly enhanced nonlinear response in composite systems, when coupled 
with a fast response time, could lead to applications in the design of switching devices 
and optical communcation systems. An excellent review on the nonlinear optical properties 
of random media has been given by Shalaev [10]. An up-to-date account on problems 
related to computational electromagnetics in dielectric heterostructures has recently been 
given by Brosseau and Beroual [11]. The current status on research on the properties of 
nanostructured random media is best reflected in the recent compilation by Shalaev [12]. 

A number of theories has been developed for weakly nonlinear composites [10,13-19], 
with the focus mainly on third-order nonlinear susceptibility such as the Kerr effect. In pre- 
vious works [20,21], we derived general expressions for nonlinear susceptibility for the second 
harmonic generation (SHG) and third harmonic generation (THG) in random composites. 
The authors also developed simple approximation for the dilute concentration limit and 
effective-medium type approximation for arbitrary concentrations. With recent advance- 
ments in the fabrication of coated particles and with possible enhanced local field effect 
in mind, we study SHG in a dilute suspension of coated spheres. Specifically, we assume 
that only the coating material has a non-vanishing SHG coefficient. The core and the host 
medium are linear in nature. We derive an expression for the effective SHG coefficient per 
unit volume of nonlinear material in the system. It is found that with modest choice of ma- 
terial parameters, the SHG coefficient can be greatly enhanced at certain frequencies. We 
also study the two-dimensional (2D) version of the problem, i.e., the SHG coefficient in a 2D 
dilute suspension of coated discs in a linear host, and similar enhancement is found. When 
the coated material has a non- vanishing THG coefficient, we give the resulting expression 
for the effective THG coefficient per unit volume of nonlinear material in both 3D and 2D 
random composites. 
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The plan of the paper is as follows. In Sec. II, we derive the effective SHG coefficient in a 
dilute suspension of coated spheres in a linear host. We illustrate the possible enhancement 
effect by considering a model system of particles with a metallic core and coated by a 
nonlinear material. In Sec. Ill, we give the corresponding results for a 2D dilute suspension 
of coated discs. Results for the effective THG coefficient in both 3D and 2D cases are given 
in Sec. IV, together with a summary. 



II. FORMALISM 

We consider the effective susceptibility for second-harmonic generation (SHG) of a dilute 
suspension of coated particles embedded in a linear host in both three dimension (3D) and 
two dimension (2D). In 3D, the system consists of coated spheres; while in 2D, the system 
consists of coated cylinders with aligned axes (or coated discs). The spheres (cylinders) are 
supposed to have inner radius r± and outer radius r 2 . The core medium is supposed to have 
(possibly complex) dielectric constant ei, the coating has dielectric constant 62, and the host 
has dielectric constant €3. In addition, the coating, but neither the host nor the core, has 
a nonlinear susceptibility for second-harmonic generation. This susceptibility is, in general, 
a tensor of rank three [22] and is denoted by dij k {—2oj\uj,uj) or dij k {uj,uj), with each of the 
subscripts running over the three Cartesian indices. Note that the elements of d may be 
complex. The total volume fraction of coating material in the composite is /, and / << 1 
is assumed. 

Hui and Stroud [20] have derived a general expression for the effective nonlinear SHG 
susceptibility dfj k (— 2u>; u, u) of a composite, under the assumption that the nonlinearity is 
weak. When applied to the system considered here, the result may be written in the form 

d ijk = f d e m n(Ku(2uj)K mj (Lj)K nk (uj)}. (I) 

Here 

K nk {x,u) = — — - (2) 

is the possible local field enhancement factor giving the n-th component of the local electric 
field at position x and frequency uj when the applied field E is in the fc-th direction at 
frequency uo. (• • •) denotes an average over the volume of the nonlinear shell and the average 
is to be calculated in the linear limit. We use the convention that repeated indices are 
summed over. 

Equations (1) and (2), as well as all the results below, are also obtained in the quasistatic 
approximation, according to which the local electric field may be written as the negative 
gradient of a scalar potential. This approximation is most accurate when the particle radii 
are small compared to both the wavelength of light in the medium and the skin depth of 
metallic particles. The quasistatic approximation neglects such effects as radiative losses, 
which may become important at shorter wavelengths; it also omits diffuse scattering, as 
discussed further in the last section of the paper. 

We consider the 3D case of coated spheres in a linear host. In the dilute concentration 
limit where the interaction between particles can be neglected, the electric field in the coated 
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shell can be calculated using standard electrostatics [23]. The electric field in the spherical 
shell E s , where the subscribe denotes the spherical shell with the coating medium, in the 
presence of an external field E z is related to the electric potential via 

E s = -V$ s (3) 

with 

$,(r,0) = -£rcos0 + £^, (4) 



where 



Es = r s E (5) 
p s = X s rlT s E . (6) 



Here T s and \ s are given by 



T = ^ ( 7 ) 

s e 2 + 2e 3 + 2(e 2 - e 3 )X 8 /fi 8 V ; 

A, = (8) 
where /i s is related to the ratio of the outer to inner radii 

lis = (r 2 / ri f. (9) 
The electric field in the coating medium can then be explicitly evaluated with the result 

_ _ , 3(xzx + yzy + z 2 z) — r 2 z , . 

E s = S s z + ^ "-JL. 1 ps? (10) 

where r is the distance from the center of the sphere to the point under consideration inside 
the spherical shell. Note that p s has dimensions of a dipole moment, and it is the dipole 
moment of the coated particle in the linear limit. Comparing with results of a single dielectric 
sphere in a host, the coated sphere behaves as a sphere with an effective dielectric constant 
e s of the form 



with \ s and ji s given by Eqs.(8) and (9). It can be immediately seen from Eq.(ll) that the 
surface plasmon resonance of a coated particle can be tuned [5,7] by the material parameters 
\ s and the geometrical parameter fi s of the particle. 

To proceed and for the purpose of illustrating the possible enhancement in SHG, we will 
make the assumption that of the possible components of d in the shell, only dm are nonzero 
and that these are all equal, i.e., dm = d{—2oj]uj,uj) = d w ^ w (i = x, y, or z). Depending on 
the symmetry of the nonlinear materials concerned, d^\. would have vanishing components 
for certain combinations of the indices. For example, one may also have d^ ^ only for 
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i ^ j Our calculations can also be carried out for the latter case. Using Eq.(l) in the 
dilute concentration limit gives 



= f d u,,u, (K zi (2u)Kl(u)) shell , 



(12) 



i=x,y,z 



where we have denoted the surviving SHG coefficient d e zzz by for simplicity. The 

brackets (• • -) s heii denote the average of the enclosed quantity over the volume of the spherical 
shell, and / is the volume fraction of coating material in the system. 

It is easily shown that only the i — z term in the above sum contributes, when a 
suitable angular average is carried out over the volume of the shell. In order to evaluate 
that remaining term in Eq.(12), we need to calculate the average over the shell volume of 
the quantity 



3z 2 - r 2 
S s (2u) + =— p s (2a;) 



3z 2 



1 2 



(13) 



where we have indicated the frequency at which the field and the dipole moment are to be 
evaluated explicitly. Defining B s = {3z 2 — r 2 )/r 5 , the average is easily evaluated by noting 
that 
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(14) 



Here v s h e ii = (47r/3)(6 3 — a 3 ) is the volume of the spherical shell. Carrying out the average 
of expression (13), the effective SHG susceptibility is then given in terms of these quantities 
by 



d e 



1 



+S 2 



S s (2uj)p 2 s (uj) + 2S s (uj) Ps (2uj)p s (uj)} + S 3 p s (2u)p 2 (u;)) 



(15) 



Equation (15) is the main result of this section. The left hand side of Eq.(15) gives the 
effective SHG susceptibility of the composite per unit volume of the nonlinear coating ma- 
terial normalized by the SHG susceptibility of the coating material. Here p s (uj) and p s (2u) 
are the dipole moments that would be induced in a linear particle when an electric field is 
applied at frequencies u or 2uj. 

The magnitude of the right hand side of Eq.(15) thus gives the possible enhancement in 
the effective SHG coefficient due to the geometry of the particles (coated particles) and the 
host medium through the factor T s . As a model system, we consider spherical particles in 
which the inner core is a Drude metal, which has a dielectric constant of the form 



+ iu/r 



(16) 
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with uj p being the plasma frequency and r the relaxation time. We take oo p = 2.28 x 10 16 s _1 
and r = 6.9 x 10~ 15 s, numbers that correspond to bulk aluminum. Note that uj p t ~ 100—200 
are typical of metals. The nonlinear shell is assumed to have a frequency independent 
dielectric constant e 2 = 2.52, while the host medium also has a frequency independent 
dielectric constant 63 = 1.76. These numbers are typical of non-conducting materials. Figure 
1 shows the real and imaginary parts of d^, u / \fdu,J) for this model system as a function of 
frequency u)/uo p and for different ratios of outer to inner radii. We note that an enhancement 
factor of the order of 10 3 can be achieved at suitable frequencies. At certain resonance 
frequency of the coated particle, the local field in the shell region is sufficiently large to drive 
a much enhanced SHG susceptibility. This frequency corresponds to the surface plasmon 
resonance of the coated particle, the value of which depends on the core material, the coating 
material, and the host medium. The frequency can also be tuned by properly selecting the 
material parameters and/or by tuning the ratio of the outer to inner radii. We see in Fig.l 
that as the ratio increases, the resonance shifts to lower frequency. We also note that higher 
enhancement per unit volume of nonlinear coating material is achieved for thinner coating. 
Therefore, one would expect a dilute suspension of spherical particles coated with a thin 
layer of nonlinear materials will give the largest enhancement per unit volume of nonlinear 
material. 



III. 2D CASE: COATED CYLINDERS 

A two dimensional random composites on the x-y plane, for example, can be considered 
as a random dispersion of cylindrical particles in a host medium, with the axes of cylinders 
aligned in the same direction perpendicular to the plane. A cross section of the system 
consists of circular discs embedded in the host. A calculation of the electric field E c in the 
coated shell of a cylinder in the presence of an applied field E x can be carried out in a way 
similar to the 3D case using standard electrostatics. Instead of Eq.(10) in 3D, we have for 
the 2D case 

_ _ ^ 2(x 2 x + xyy) — r 2 x , , 

E c = S c x + ^ p Pc , (17) 

where r is the distance from the center of a disc (or cylinder) to the position under consid- 
eration inside the coating shell. Here, 

E C = T C E , (18) 
Pc = \ c rjT c E , (19) 

Fc = I J"" U / ' ( 2 °) 
£2 + £3 + (e 2 - e 3 )A c //i c 

+ e 2 

He = (-) 2 , (22) 

n 

and the subscript "c" is used to indicate that the case of coated cylinders or discs are being 
considered. Comparing with results of a single dielectric circular disc in a host, a coated 
disc behaves as a disc with an effective dielectric constant e c of the form 
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_ _ 1 + \c/Hc 

e c — ~ 7—, — £2, 

1 - Xc/fJ-c 



(23) 



with A c and /jl c given by Eqs.(21) and (22). 

To obtain an expression for 2D similar to Eq.(15), we make the same assumption for the 
SHG coefficient of the coating material. This leads to an expression similar to Eq.(12) for 
the effective SHG suscepbility, except that we need to evaluate the average over the shell 
area of the quantity 





S c {2u) + Pc (2uj) 



£ M + — i — Pc(^) 



(24) 



instead of the expression (13) in 3D. Defining B c = (2x 2 
evaluated by noting that 



r 2 )/r 4 , the average is easily 
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(25) 



Here a she u is the area of the circular shell (ring) of coating. 
The effective SHG susceptibility is then found to be 



[£ c (2u)£ 2 (u) + S 2 {£ c {2u)p 2 c {u) + 2£ c {u)p c {2u)p 2 c {u))\. 



(26) 



Equation (26) is the result in 2D, analogous to Eq.(15) in 3D. It gives the effective SHG 
suscepbitility per unit area of nonlinear coating medium normalized to the SHG susceptibility 
of the coating medium in a 2D dilute suspension of coated discs in a linear host. 

Figure 2 shows w / (fd UtU) ) as a function of frequency in the 2D case, for the same model 
system considered in the previous section. The features are similar to those in the 3D case. 
The enhancement factor tends to be larger in 2D than in 3D for the same set of model 
parameters. For discs coated with thin layer of nonlinear material, the SHG response per 
unit area of nonlinear material can be driven by the local field to achieve an enhancement 
of 10 3 , for modest choice of material parameters. 



IV. DISCUSSION 

It is straightforward to generalize the formalism to third harmonic generation (THG) 
[21]. For coating nonlinear material in which the leading nonlinear response is THG with 
the THG coefficient denoted by Xw,u,wi we can again use the local electric field at the shell 
region to drive an enhanced THG response per unit volume of nonlinear material in a dilute 
suspension. The derivation is similar to that in Sec. II, although somewhat tedious. The 
result is 
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[£ s (3u})£s(u) + 3S 2 £.s(uj)Ps(uj)(£ s (uj)p s (3uj) + £ s (3uj)p s (uj)) + 



£b(3w).E§(c<;) 

+S 3 ^H(3£ s (u;)p s (3u;) + £ s (3u;)p») + ^(3^)^))], (27) 

where 

^ = (28) 

for a dilute suspension of coated spherical particles. 

In the 2D case of coated discs of nonlinear THG coating material, the result is 



\£ c {3u)£l{u) + 3S 2 £Mp c {uj){£ c {uj)p c {3u) + £ c {3u)p c {u)) + 



+S 4Pc (3uj)p 3 c (uj)}, (29) 



where 



as/ie«« r x r 2 8 



To illustrate the enhancement in THG susceptibility, we use the same choice of material 
paramters as our model system and calculate the enhancement factor corresponding to 
Eq.(27) for spherical particles. Results are shown in Fig. 3. Since the THG is driven by 
higher power of the local field, the enhancement factor is found to be of the order of 10 5 at 
resonance frequency, for our choice of parameters. 

Note that, when we calculate the effective SHG and THG susceptibilities in the present 
work, we do so by treating the medium as effectively homogeneous, with effective properties. 
Thus, our calculation would lead to a collimated beam of harmonic emission when this effec- 
tive medium is used in an SHG experiment. Besides this collimated emission, some diffuse 
scattering at the second-harmonic frequency is also expected, because of the inhomogeneity 
of the medium. Indeed, such diffuse SHG scattering has been observed in a recent experi- 
ment, carried out near the percolation threshold of a random metal-insulator composite [24], 
in addition to the collimated beam described here. This diffusive scattering should have an 
intensity varying as a power of the particle radius. The present approach would need to be 
extended to include such diffuse scattering. 

Finally, we briefly comment on the validity of the low-concentration approximation used 
to obtain eqs. (15), (16), (27) and (29). These equations are obtained by calculating the 
fields and displacement vectors as if the coated particles feel a uniform applied field. In ac- 
tuality, the local field should include the dipole fields produced by the other particles. These 
fields may start to become significant when the dipole field due to one particle, calculated 
at the position of a neighboring particle, becomes comparable to the applied electric field. 
Although this condition may be satisfied near resonance even for / <C 1, the large enhance- 
ments we find may still persist at higher /, because the interaction fields, even though they 
can be substantial, may approximately cancel out if the environment of a particle can be 
approximated as isotropic. The local field factors at higher concentrations can be estimated 
using such methods as the Maxwell-Garnett or effective-medium approximations. 
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Our results for SHG (Eqs.(15) and (26)) and THG (Eqs.(27) and (29)) are applicable 
to other possible structures of coated particles, as long as the composite is in the dilute 
concentration limit and SHG originates from the coating material. One could use a dielectric 
core and a metallic shell with a non- vanishing SHG coefficient [25] . In general, careful choices 
can be made so that the local electric field in the coating region is enhanced. This enhanced 
local field can then be used to drive the SHG response. In addition, it is possible to tune 
the local field by tuning the core size and the shell thickness. For composites of higher 
concentrations, the Maxwell-Garnett approximation and the effective-medium appoximation 
can be invoked in estimating the local field factors. In higher concentrations, percolation 
effects can also lead to possible enhancement in SHG and other nonlinear response. 
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FIGURES 



FIG. 1. The real and imaginary parts of the effective SHG coefficient d^, u /{fdw,u) P er vna.it 
volume of nonlinear coating material as a function of frequency to/ujp for a dilute suspension of 
particles with metallic core coated with a nonlinear material. Results for three different ratios of 
the outer to inner radii fi s of the particles are shown. 

FIG. 2. The real and imaginary parts of the effective SHG coefficient d% ^/{fd^^) per unit 
area of nonlinear coating material as a function of frequency lu/lo p for a dilute suspension of discs 
with metallic core coated with a nonlinear material in a 2D system. Results for three different 
ratios of the outer to inner radii fj, c of the discs are shown. 

FIG. 3. The real and imaginary parts of the effective THG coefficient xZ,uj,uj/ ifXu,u,u) P er un it 
volume of nonlinear coating material as a function of frequency uj/lo p for a dilute suspension of 
particles with metallic core coated with a nonlinear material. Results for three different ratios of 
the outer to inner radii fi s of the particles are shown. 
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nonlinear shell and metal core ( 2D ) 
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